Density dependence of the "symmetry energy" in the lattice gas model 
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Isoscaling behavior of the statistical emission fragments from the equilibrated sources with Z = 30 
and N = 30, 33, 36 and 39, resepectively, is investigated in the framework of isospin dependent lattice 
gas model. The dependences of isoscaling parameters a on source isospin asymmetry, temperature 
and freeze-out density are studied and the "symmetry energy" is deduced from isoscaling parameters. 
Results show that "symmetry energy" C sy m is insensitive to the change of temperature but follows 
the power-law dependence on the freeze-out density p. The later gives C sym = 30(p/po) ' 62 if 
the suitable asymmetric nucleon-nucleon potential is taken. The effect of strength of asymmetry 
of nucleon-nucleon interaction potential on the density dependence of the "symmetry energy" is 
dicussed. 

PACS numbers: 24.10.Pa, 25.70.-z, 05.70.Ce 
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I. INTRODUCTION 

One of major goals of nuclear physics is to understand 
nature of the nuclear matter, such as the nuclear equa- 
tion of states (EOS) and interaction of their basic bull- 
ing blocks, namely neutrons and protons. Especially, the 
propoerties of the asymmetric nuclear matter is of great 
interesting in recent years since it has important impli- 
cations for astrophysics and nuclear physics studies. The 
availability of beams with large neutron-to-proton ratio, 
N/Z, provides the opportunity to explore the symmetry 
energy term in nuclear equation of state in very asym- 
metric nuclear systems. In such reactions, isospin degree 
of freedom has a prominent role and can serve as a valu- 
able probe of the symmetry energy term of the nuclear 
equation of state. Recently the isoscaling behavior ex- 
tracted from the fragment isotopic composition between 
two reactions with similar entrance channel has been pro- 
posed as a very powerful tool to determine the symmetry 
energy term of EOS This scaling law relates ratios of 
isotope yields measured in two different nuclear reactions, 
1 and 2, R 21 (N,Z) = Y 2 (N, Z)/Y 1 (N, Z), and such ra- 
tios are shown to obey an exponential dependence on the 
neutron number N or proton number Z of the isotopes 
or isotones characterized by three parameters a, (3 and 
C0: 



R 21 {N, Z) 



Y 2 (N, Z) 
Yi(N, Z) 



Cexp{aN + PZ), (1) 



where C is an overall normalization constant. So far, 
this universal law expression for R 2 \ has been extensively 
explored experimentally and theoretically, under diverse 
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approximations, to primary yields produced by disassem- 
bling infinite equilibrated systems in microcanonical and 
grand canonical ensembles as well in canonical ensembles. 
It has been observed in statistical models 0, [H, 0] , lattice 
gas model [5| , percolation model [6| , statistical abrasion- 
ablation model 0, 3, fission model 0, as well as the 
molecular dynamics type models [Tol . fill ] . Furthermore, 
under these approximations, the isoscaling parameters a 
and (3 have been found to be related to the symmetry 
term of the nuclear binding energy [l], 0] , to the level of 
isospin equilibration (T^ |. and to the values of transport 
coefficients [l3j . 

The paper presents results of the isoscaling analy- 
sis for the light fragments from thermal sources which 
are produced by the lattice gas model. The tempera- 
ture and density dependences of of the symmetry energy 
are extracted, and the effect of interaction potential be- 
tween nucleon-nucleon on the symmetry energy is also 
discussed. 

The article is organized as follows. Section [IT] makes a 
simple introduction on the lattice gas model. In Section 
IIII1 isoscaling phenomenon and its dependence on tem- 
perature and density are presented, while the tempera- 
ture and density dependence of the dedced "symmetry 
energy" from the isoscaling parameter has been deter- 
mined. The effect of asymmertic nucleon-nucleon poten- 
tial on the "symmetry energy" is also discussed. Finally 
a summary is given in Sec. IV. 



II. MODEL OVERVIEW 

The lattice gas model was developed to describe the 
liquid-gas phase transition for atomic system by Lee and 
Yang [14]]. The same model has already been applied to 
nuclear physics for isospin asymmetrical systems jlH [l6| . 
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In the lattice gas model, A (= N + Z) nucleons with an 
occupation number S{ which is defined Si = 1 (-1) for a 
proton (neutron) or s, = for a vacancy, are placed on 
the L sites of lattice. Nucleons in the nearest neighboring 
sites interact with an energy e SiSj . The hamiltonian is 
written as 



E 



A 

£ 

i=l 



P 2 
2m 



Kj 



. SiS j 



(2) 



In order to investigate the symmetrical term of nuclear 
potential in this model, we can use sets of parameters: 
one is an attractive potential constant e SiSj between neu- 
tron and proton but no interaction between like nucleon, 
i.e. proton and proton or neutron and neutron, namely 
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£nn = £pp = O.MeV, 

e pn = -5.33MeV. (3) 

This kind interaction results in an asymmetrical poten- 
tial among different kind of nucleon pairs, hence it is an 
isospin dependent potential. In previous studies, e pn is 
generally chosen to be -5.33 MeV [15j, [l6[. Later we will 
use different {e pn ,e nn ) parameters to compare the effect 
of this asymmetrical potential on the isocaling parame- 
ters and "symmetry energy" . 

In the LGM simulation, a three-dimension cubic lattice 
with L sites is used. The freeze-out density of disassem- 
bling system is assumed to be pf = j;po, where po is the 
normal nuclear density. The disassembly of the system is 
to be calculated at pf , beyond which nucleons are too far 
apart to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum 
by Monte Carlo samplings of Maxwell-Boltzmann distri- 
bution. Once this is done the LGM immediately gives 
the cluster distribution using the rule that two nucle- 
ons are part of the same cluster if their relative kinetic 
energy is insufficient to overcome the attractive bond, 
i.e. — e SiSj SiSj < 0. This method is similar to 

the Coniglio-Klein's prescription [T3| in condensed mat- 
ter physics. Then the fragment isotopic composition can 
be acquired. 

In the lattice gas model the ground state energy (bind- 
ing energy) can be defined as 

Ebind — £pnNp n -\- C nn N nn , (4) 

where N pn or N nn is the number of the bonds between 
unlike nucleons or like nucleons in the ground state which 
corresponds to a cold nucleus at zero temperature and 
normal nuclear density in our present classical model . 
Practically, N np or N nn is determined by the geometry 
and is equal to the maximum bond number of unlike 
nucleons or like nucleons of the systems. 



FIG. 1: (Color online) The ratio of the isotopes for Z < 9 
from two different sources for A = 66 and 63 with Z — 30 
at T = 3 MeV and pf = 0.48po as a function of the neutron 
number N. 



III. RESULTS AND DISCUSSIONS 

A. Isoscaling behavior 

To make a systematic study of source parameters 
which might influence the isoscaling behavior, several 
pairs of equilibrated sources are considered at various 
initial temperature T and freeze-out density pf in our 
present work. The equilibrated source pairs are chosen 
with different isospin asymmetry N/Z with the fixed Z 
= 30 but A = 60, 63, 66 and 69, respectively, which cor- 
responds to N/Z = 1, 1.1, 1.2 and 1.3, respectively. We 
adopt the widely used convention to denote with the in- 
dex "2" the more neutron-rich system and with the index 
"1" the more neutron-poor system. In this situation the 
value of a is always positive because more neutron-rich 
clusters will be produced by the neutron-richer source 
and the value of (3 is always negative. The yield ratios 
Rzx{N,Z) are calculated and the corresponding isoscal- 
ing behaviors are investigated over all possible decayed 
fragments. The cubic lattices of 5 3 , 6 3 , 7 3 , 8 3 and 9 3 
are used for all systems, which gives the freeze-out den- 
sity p f = 0.48 po, 0.28 po, 0.18 p , 0.12 p and 0.08 p , 
respectively, for A = 60 system. There are slight differ- 
ent freeze-out densities for A = 63, 66 and 69 systems. 
100000 events were simulated for each T-pf combination 
of each source which ensures good statistics for results. 
In the present study, we will focus on the light fragment 
isotopic distribution to explore the isoscaling behavior 
and symmetry energy. The advantage of the LGM is 
that the source density and temperature parameters are 
well defined. 

Fig. [T] shows an example of the isoscaling behavior. It 
plots the i?2i as a function of fragment neutron number 
N from which the slope parameter (isoscaling parameter) 
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FIG. 2: (Color online) (a) Temperature and density depen- 
dences of isoscaling parameters a extracted from source pair 
with A = 63 and 60. (b) Same as (a) but for the deduced 
"symmetry energy". 

can be extracted. Average isoscaling parameter a which 
is calculated over the range Z < 9 is used to discuss the 
dependence of a on the properties of emission source, 
such as temperature, source size and source asymmetry 
of the isospin. In the following calculations, isoscaling 
parameters a refer to the calculated average values. 

B. Temperature and density dependence of 
symmetry energy 

To explore the origin of isoscaling behavior and the de- 
pendence of isoscaling parameters a on the isospin com- 
position, we performed calculations on source systems 
with different densities and source asymmetry (N/Z) val- 
ues. The temperature and density dependences of a from 
the source pair (N = 33, Z = 30) and (N = 30, Z = 30) 
are shown in Fig. Ufa). Here the scattering points are 
the calculation results for the LGM with asymmetrical 
nucleon-nucleon potential, i.e. Eq.(3). Note that the sim- 
ilar behavior for other source pairs was obtained but not 
shown here. From the figure, a decreases as the temper- 
ature arises, which indicates that the isospin dependence 
of the fragment yields becomes weaker with the increas- 
ing of temperature. But an increasing trend of a with 



FIG. 3: (Color online) (a) a as a function of freeze-out density 
for various source pairs with Z — 30 which is illustrated in 
the figure at temperature T = 3 MeV. (b) The deduced "sym- 
metry energy" as a function of freeze-out density. The line 
represents the fits with C cy m = 24.1(p/po)°' 59 - Symbols in fig- 
ure correspond to YX=63/YX=60 (solid square), YX=66/Ya=63 
(open circles) and Ya=69/Ya,=66 (solid triangles). 

density is also clearly seen from the figure, which shows 
that lower density system has lower chemical potential 
that makes smaller yield probable. 

It has been shown in the framework of the grand- 
canonical limit, a has the form, 

a = ^pA[(|) 2 ] = ^p[(§)? - (§)!], (5) 

where C sym is the deduced "symmetry energy" coeffi- 
cient (MeV), (4)f = i 2 means the square of Z over A for 
system 1 and 2. T is the temperature of the system 
in MeV. Above equation has been proved to be a good 
approximation by many calculations and experimental 
data, and are generally adopted to constrain the symme- 
try energy coefficient C sym in experiments. This behav- 
ior is attributed to the difference of isospin asymmetry 
between two reaction systems in similar nuclear temper- 
ature. Since the symmetry energy has an important role 
on nuclear structure of neutron-rich or neutron-deficient 
rare isotopes, studies on the isoscaling behavior can be 
used to probe the isospin dependent nuclear equation of 
state. Fig. (Hb) shows that the deduced "symmetry en- 
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ergy" C sym almost keep flat regardless the change of tem- 
perature. In other words, the "symmety energy" shows 
little dependence on temparture. This is consistent with 
the analysis of the symmetry energy versus the temper- 
ature within a degenerate Fermi gas model [L9J. Fur- 
thermore, the deduced "symmetry energy" displays an 
increasing trend with the freeze-out density, which will 
be discussed later. 

Fig. [3|a) shows that the a increases with the density 
for the three source pairs at T = 3 McV. If the Eq. (4) 
is satisfied in our simulation, a has a linear dependence 
on A[(^) 2 ]. Since this parameter is dependent on the 
source systems, we divide a by A[(^-) 2 ] to remove the 
system isospin and size dependences. Furthermore, since 
we know value of temperature for each simulation, we 
can use the Eq.(4) to extract "symmetric energy" pa- 
rameter. Fig. EJb) shows the "symmetric energy" pa- 
rameter Csym — aT/[4A(^) 2 ] extracted from the ratios 
of different reaction system pairs, namely Ya=6q/Ya 3 =6q, 
^=66/1^=63, and Y a - S3 /Ya=60, as a function of density. 
Obviously, the same dependence is displayed even though 
for the different source combinations. Considering that 
the generally used density-dependent symmetry energy 
coefficient has been taken, namely 

C cym = C (p/p ) J , (6) 

we can use the above expression to fit C sym (p) . The line 
in the Fig.^b) shows the fit. It gives C ~ 24 MeV and 
7 ~ 0.6. It was worthy to mention that 7 value is in a 
reasonable range, namely 0.55-0.69, which was proposed 
by recent data and models [20| . but the Co looks softer 
than the recent suggested value [2(|, namely 31-33 MeV, 
even though it is consistent with the expanding emitting 
source model calculatiopn by Tsang et al. [l|. However, 
it must be mentioned that the binding energy per nucleon 
Ebind/A from the Eq. 2] for the studied system in our 
LGM with (e pn = -5.33, e n „ = 0) is around -11.7 MeV 
which is greater than -16 MeV, the energy per nucleon in 
a uniform nuclear matter at T = at the normal density 
Pq. Later we will find the potential parameter set (e pn = 
— 7, e n n = 0) will give a more reasonable binding energy 
as well as Co. 

C. Symmetry energy and nucleon- nucleon 
symmetric potential 

In order to discuss the effect of the isospin dependent 
interaction potential on the deduced "symmetry energy" , 
we plot Fig. 2] for T — 3 MeV, by changing the symmet- 
ric potential, i.e. on the one hand we change e pp from 
0, -1, -3 to -5.33 MeV but keep e pn = -5.33 MeV con- 
stantly; on the other hand we change (e pn — -7, e pp — 
0) to comparison with the normal case, namely e pn — 
-5.33 and e pp = 0. In fact, these five (e pn , e pp ) configura- 
tions, namely from (-7,0), (-5.33,0), (-5.33,-1), (-5.33,-3) 
till (-5.33,-5.33), counterparts gradually decreased asym- 
metric degree which can characterized by the value of 




FIG. 4: (Color online) Deduced "symmetry energy" for (N 
= 33, Z = 30) and (N = 30, Z = 30) source pairs at T 
= 3 MeV as a function of density with different (epn, e pp ) 
parameters. The detailed values for e nn and e pn are indicated 
in the figure. Lines represent the fits with Csym = Co(p/po) 7 
and the corresponding detail is illustrated in the figure. 



e pn -e pp . When the symmetric potential becomes not sen- 
sitive slowly to nucleus type, that is close to homogeneity. 
It is observed that the isoscaling property is preserved 
for all situations, but the values of the deduced "sym- 
metry energy" parameters shows different sensitivity to 
the freeze-out density in different condition of (e pn , e P v). 
As derived in other models and theoretical frame [3, HJ , 
isotope yield ratio is dominantly determined by the sym- 
metry term in the binding energy for two equilibrium 
sources with comparable mass and temperature but dif- 
ferent isospin degree N/Z. It is clear that the values of 
the deduced " symmerty energy" in the isoscaling relation 
are related to the symmetry potential between nucleon 
and nucleon: when nucleon- nucleon has smaller symme- 
try potential (i.e. smaller difference between e pn and e pp ) 
of the system, the value of C sym decreases. To quantita- 
tively see the density dependence of "symmetry energy" , 
we fit the points with the Eq.(6). The parameters of Co 
and 7 are shown in the figure. Generally, both Co and 7 
increase with the asymmetry of the nucleon-nucleon po- 
tential. In other words, the "symmetry energy" param- 
eter Co tends to increase and asymmetric EOS becomes 
harder with the asymmetry of the nucleon-nucleon poten- 
tial. It is of very interesting to see Co ~ 30 MeV and 7 ~ 
0.62 for the potential parameter set (e pn — 7,e nn = 0). 
The above parameters are exactly consistent with the re- 
cently suggested form of density dependence of the sym- 
metry energy [2(| ■ This is not accidental since the LGM 
with the potential parameter set (e pn = 7, e nn = 0) gives 
a reasonable binding energy -15.4 MeV/A for the stud- 
ied systems, which is close to the energy in a uniform 
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nuclear matter at T = at the normal density po- 

IV. CONCLUSION 

In summary, the isotopic effect of the light fragments 
from the equilibrated thermal sources for the source pairs 
N = 30, 33, 36 and 39 with constant Z = 30 has been 
successfully investigated by the lattice gas model model. 
The model constrains the source temperature and frag- 
menting density which is below saturate nuclear density 
Po, and the fragmentation starts from an equilibrated 
system which is separated from any dynamical process. 
In our work we change input parameters of tempera- 
ture and freeze-out density to survey the evolution of 
isospin degree of freedom on the isoscaling behavior and 
the symmetry energy. Isoscaling phenomenon is observed 
for the emitted light fragments and their dependences on 
source isospin asymmetry, temperature and density have 
been systematically investigated. Isoscaling parameters 
a decrease with the increasing of temperature as well 
as the decreasing of density. The deduced "symmetry 
energy" C sym was extracted from the isoscaling parame- 



ters, which shows little dependence on temperature but 
follows the power-law dependence on the freeze-out den- 
sity, namely C sym — Co(p/po) 7 - It is shown that Co ~ 30 
MeV and 7 ~ 0.62 in case that the asymmetry nucleon- 
nucleon potenatial (e pn = -7, e pp = 0) is used, this result 
is consistent with the recently suggested values from data 
and models. Finally the effect of the asymmetry nucleon- 
nucleon potenatial on the deduced "symmetry energy" 
are studied by adjusting (e pn , e pp ) parameters, it shows 
that Co becomes larger and the density dependence of 
symmetry energy becomes stiffer with the increasing of 
asymmetric part of the nucleon-nucleon potential. There- 
fore it is a powerful tool to use the isoscaling method to 
extract the asymmetry EOS from the confrontation of 
data with the models. 
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